Cooper pairing and BCS-BEC evolution in mixed-dimensional Fermi gases 
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Similar to what has recently been achieved with Bose-Bose mixtures [1], mixed-dimensional Fermi-Fermi 
mixtures can be created by applying a species-selective one-dimensional optical lattice to a two-species Fermi 
gas (cr = {f, 4-}), such a way that both species are confined to quasi-two-dimensional geometries determined 
by their hoppings along the lattice direction. We investigate the ground state phase diagram of superfluidity for 
such mixtures in the BCS-BEC evolution, and find normal, gapped superfluid, gapless superfluid, and phase 
separated regions. In particular, we find a stable gapless superfluid phase where the unpaired f and 4- fermions 
coexist with the paired (or superfluid) ones in different momentum space regions. This phase is in some ways 
similar to the Sarma state found in mixtures with unequal densities, but in our case, the gapless superfluid phase 
is unpolarized and most importantly it is stable against phase separation. 

PACS numbers: 03.75.Ss, 03.75.-b 



I. INTRODUCTION 

Atomic Fermi gases have emerged as unique testing ground 
for many theories of exotic matter in nature, allowing for the 
creation of complex yet very controllable many-body quan- 
tum systems [2], where for instance observation of the BCS- 
BEC crossover has so far been the most important achieve- 
ment in this field. Following this huge success with single- 
species fermion mixtures, there has been increasing exper- 
imental interest in studying two-species Fermi-Fermi mix- 
tures JH-d] . In particular, 6 Li- 40 K mixtures have recently been 
trapped and interspecies Feshbach resonances have been iden- 
tified, opening a new frontier in ultracold atom research to 
study exotic many-body phenomena, one of which is the pos- 
sibility of studying fermion pairing in mixed dimensions [9- 

Eh. 

Mixed-dimensional atomic systems, in which two types of 
particles live in different dimensions, can be created with 
two-species Fermi-Fermi, Bose-Fermi, and Bose-Bose mix- 
tures by using species-selective optical lattice potentials. This 
has recently been achieved with Bose-Bose mixtures |Q]] by 
applying a one-dimensional optical lattice to the 41 K- 87 Rb 
mixture, where only 41 K atoms feel the lattice potential, and 
they are confined to a quasi-two-dimensional geometry, while 
having negligible effect on 87 Rb atoms, that is leaving 87 Rb 
atoms three dimensional. Motivated by this experimental 
work, here we analyze Cooper pairing in mixed-dimensional 
Fermi gases. We consider both single-species and two-species 
fermion mixtures, and analyze the ground state phase dia- 
grams in the BCS-BEC evolution, which involves normal, 
gapped superfluid, gapless superfluid, and phase separated re- 
gions. In particular, the gapless superfluid phase, where the 
unpaired j" and J, fermions coexist with the paired (or super- 
fluid) ones in different momentum space regions, is in some 



ways similar to the Sarma state found in mixtures with un- 
equal densities lfl2ll . but in our case, the gapless superfluid 
phase is unpolarized and most importantly it is stable against 
phase separation. In this way, our gapl ess superfluid phase 
is very similar to those of Refs. 11131 U4I1 . which are recently 
proposed for ultracold atomic systems in other contexts. 

The rest of the manuscript is organized as follows. In 
Sec. lU after introducing the Hamiltonian in Sec. Ill Al the 
corresponding saddle-point self-consistency equations are de- 
rived in Sec. Ill Bl and their noninteracting limit is discussed in 
Sec. Ill CI We numerically solve these equations in the BCS- 
BEC evolution in Sec. [TTTJ where we investigate the normal- 
superfluid transition in Sec. IIII Al the topological gapless- 
gapped superfluidity transition in Sec. IIII B I and the ground 
state phase diagrams in Sec. IIII CI A brief summary of our 
conclusions is given in Sec. [IV] We also include three ap- 
pendices, where the self-consistency equations are further dis- 
cussed in Appendix A, boundary equation for the normal- 
superfluid transition is derived in Appendix B, and the molec- 
ular BEC limit is investigated in Appendix C. 



* Current address: 47'nci Muht. Bl. K.ligi PK:75 Selcuklu, Konya, Turkey. 



II. MIXED-DIMENSIONAL FERMI GASES 



In this work, we analyze Cooper pairing in mixed- 
dimensional Fermi gases, which seems to be a very promising 
way to create superfluidity with mismatched Fermi surfaces, 
and the physics involved is in some ways similar to that of the 
unequal density mixtures Il5l - [l8ll . We consider only uniform 
(homogenous) mixtures, but emphasize that the finite-size ef- 
fects due to the confining trapping potentials (which are al- 
ways present in atomic systems) can be taken into account 
using the local-density approximation (as a first approxima- 
tion). 
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A. Hamiltonian 

To describe such mixed-dimensional Fermi gases in a 
species-selective one-dimensional optical lattice (say in the 
z direction), we start with the real-space Hamiltonian (h = 

k B = 1) 



H = E / d " r <o 



2m 



(1) 



where the pseudo-spin a labels both the type and hyperfine 
states of atoms represented by the creation operator tpl a , 
and m a is the mass and [i a is the chemical potential. Here, 

r = (r x ,r y , r z ) is the position with rj_ — ^Jr 2 + r 2 , Vq L = 

Vb, CT sin 2 (7rr z /(i z ) is the optical lattice potential, d z is the lat- 
tice spacing, and g > is the strength of the attractive inter- 
action (zero-range and isotropic) between j" and J, fermions. 

In order to achieve the momentum-space Hamiltonian, we 
first expand the creation and annihilation field operators in 
the orthonormal and complete basis set of Wannier func- 
tions W(r z — r l z ) of the lowest-energy states of the op- 
tical potential near their minima (single-band approxima- 
tion), and then take the Fourier transform of the site opera- 
tors. This corresponds to the following substitution: V^r.cr = 
(1/VM) E i)k a K aW a (r z - ri)e- i ^ T ^ +k ' r i\ where k = 
(k x , k y ,k z ) is the momentum with k± 



'k 2 + ky, i labels 

the lattice sites, and M is the total number of them. Keep- 
ing only the tunneling between nearest-neighbor sites and on- 
site interactions, and using the orthonormality of the Wannier 
functions, the resultant Hamiltonian can be written as 



H = 



k,<r 



t t 



9 



&k',q> 



(2) 



k,k',q 



where & fc q = a k+q / 2 t a -k+q/2 J, creates fermion pairs 
with center of mass momentum q and relative mo- 
mentum 2k. Here, £k,<r = ek,er — AV, where 
£k.cr = k±/(2m cr ) + 2t a [1 — cos(k z d z )] is the single- 
particle energy dispersion, and t a — J dr z W*(r z — 
4) [-d 2 /(2m a dr 2 z ) + V§ L {r z )] W(r z - r{) is the tunnel- 
ing amplitude between any nearest-neighbor sites i and j. 

Following the usual treatment, strength of the attractive 
interaction can be written in terms of an "effective" s- 
wave scattering length a e ff as 1/g = —m + V/(4:ira e ff) + 
^ k l/(2e kj+ ), where m± = 2m^m^j [m^ ± m^), V is the 
volume of the system, and ek,± = (fik.t ± £k,4,)/2. Note that 
m + is twice the reduced mass of the f and \. fermions, and 
that the equal mass case corresponds to |to_| — > oo. Here, 
and throughout, the momentum space sums are evaluated as 
£ k = [V/{2nf] J d*k ee [V/(2n 2 )} fj d ~~ dk z f Q °° k ± dk ± , 
since the system at hand has a cylindrical symmetry around 
the k z -axis, and a translational symmetry along the z direc- 
tion, so that the k z integrals are limited to the first Brillouin 
zone, i.e. —ir/d z < k z < ir/d z . The resultant integrands are 



also even functions of k z , and hence we integrate over half of 
the Brillouin zone and multiply them by two. 

In this manuscript, for its simplicity, we set tj, = 
1 / {2m] r d 2 ) but allow for the f fermions to have a different ef- 
fective mass along the lattice direction through a tight-binding 
dispersion determined by tf, so that 



£k,t 



2m-f 
k\ 



2m; m^d 2 z 



2t t [1 - cos(fc z <4)] . ( 3 ) 



[1 - cos{k z d z )} . (4) 



Notice that when kp z ±d z <C ir (i.e. the low-filling limit 
for the I species), where kp z ^ is the Fermi momentum of 
I fermions in the k z direction, since Eq. (O can be approxi- 
mated as ek,j, ~ k 2 / (2to^), optical lattice has negligibe effect 
on 4- fermions in this limit. Next, we analyze the phase dia- 
gram of the Hamiltonian given in Eq. (O with the dispersions 
given by Eqs. (0 and |@), within the saddle-point (mean-field) 
approximation. 



B. Self-consistency equations 

At low temperatures (T w 0), the saddle-point self- 
consistency (order parameter and number) equations are suf- 
ficient to describe the BCS-BEC evolution of superfluid- 
ity lfl9l l20ll . For the Hamiltonian given in Eq. (O, the saddle- 
point action is So = fio /T, where 

r> = ^- + ]T (&,+ - E K+ ) + TJ2 ln[/(-£ k>s )](5) 

" k k,s 

is the saddle-point thermodynamic potential. Here, f(x) = 
l/[exp(a;/r) + 1] is the Fermi function, £k,s = (£ k + + 
lAp) 1 / 2 + 7s^k,- is the quasiparticle energy when 71 = 1 
or the negative of the quasihole energy when 72 = — 1, and 
-E k ,± = (£k,i ± Ek t 2)/2. In addition, A is the order parame- 
ter and £k,± = £ k,± — M±j where n± = (^±^)/2. Note that 
the symmetry between quasiparticles and quasiholes is broken 
when £ k ._ 7^ 0. 

The saddle -point condition <5S!o / 5A* = leads to an equa- 
tion for the order parameter, 



4ira 



E 



1 - /(£k,i) - f{E K 2) 1 



2E ] 



k,+ 



, (6) 



where, as usual, g is eliminated in favor of the effective s- 
wave scattering length a e j / via the relation given above lfl9l 
l20ll . The order parameter equation has to be solved self- 
consistently with the number equations. At the saddle point, 
the relation = — dH,o/d/j, a leads to 

JV t = E [l u k| 2 /(^ka) + M 2 /(-£k, 2 )] , (7) 

k 

Ni = [l«k| 2 /(^k, 2 ) + bk| 2 /(-^k,i)] , (8) 
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where |u k | 2 = (1 + £ k ,+/£ k ,+)/2 and \v k \ 2 = (1 - 
£k.+/£'k,+)/2 are the usual coherence factors. At T = 0, 
Eqs. (O, (0 and ([8]) can be simplified considerably as shown 
in Appendix A. 

In order to analyze the phase diagram at T = 0, we solve 
the saddle-point self-consistency equations and check the sta- 
bility of these solutions for the uniform superfiuid phase us- 
ing the compressibility (or the curvature) criterion I2H - I23I1 . 
This says that the compressibility matrix k(T) with elements 
k<x,ct'(T) = —d 2 no/(d/j,o-dfi c7 ' ) needs to be positive definite, 
and it is related (identical) to the condition that the curvature 



d 2 n 

~dA 2 



IAI 



E 



0.5- f(E Ks ) , f'(E ktS ) 



E 3 



E 2 



(9) 



of the saddle-point thermodynamic potential flo with respect 
to the saddle-point parameter A needs to be positive. Here, 
f'(x) = df(x)/dx. When at least one of the eigenvalues of 
k(T), or the curvature d 2 £lo/dA 2 is negative, the uniform 
saddle-point solution does not correspond to a minimum of 
£lo, and a nonuniform superfiuid phase, e.g. a phase separa- 
tion, is favored. 



C. Noninteracting Limit 

Before presenting our numerical results, let's first analyze 
the noninteracting j->0 limit. In this limit, since 1 / a e f f — > 
—00, the order parameter vanishes A — > 0, and at T = 
Eqs. (|7} and © can be written as N a = ^ k 9(— £k,o-)7 where 
8(x) is the heaviside step function. Evaluating the k-space 
sums, the density n a = N a /V of the a fermions become 



m a k Fz 



-(AV-2< CT ) + 



2ir 2 ir 2 d z 

Ha = 2t a [1 - COs(k Fz ,crd z )] , 



■sm(k Fz ^d z ), (10) 
(11) 



t,k 2 d 2 



where k Fz ^ a > corresponds to the Fermi momentum of the 
a fermions in the k z direction, which is related to the chemical 
potential via Eq. (fTTT) . 

Note that, in the low-density (k Fz ^d z <C if) limit, expand- 
ing out Eqs. ( [Tot and (fTTT i to the lowest nontrivial orders (third 
and second orders, respectively) in k Fz ^d z , we obtain n 
m^kpz^i ~ h k F Z ,i d z/3)/( 27T2 ) and ^ 
Combining these two expressions, and setting k Fz ^ = kp.i 
or ti = l/(2mid 2 ), the density of I fermions acquires the 
usual form ss k F ,/ (6ir 2 ). In addition, when — > and 
&Fz,t — 7T /d z , Eq. ( fTOb reduces to the density of j" fermions 
in each (isolated) two-dimensional planes along the k z direc- 
tion, N-\/(MA) — m-fkp ^ j '(47r) , which is of the usual form, 
where we used V — AL Z with A being the area of the system 
in the (x, y) plane and L z is the system size in the z direction. 
Here, M = L z /d z is the number of two-dimensional planes, 
i.e. the number of lattice sites, in the z direction. Having 
discussed the noninteracting limit, next we analyze the BCS- 
BEC evolution. 



III. SADDLE-POINT APPROXIMATION 



In this section, we consider only equal density mixtures, 
where — nu at zero temperature. For these mixtures, first 
we analyze the amplitude of the order parameter | A|, chemical 
potential sum and the chemical potential difference \i_ as 
a function of the tunneling amplitude t-j- and effective scat- 
tering length l/(fci^a e //) for fixed values of k F ^d z . Here, 
k F ,i is an "effective" Fermi momentum for \, fermions defined 
through the three-dimensional density = k F ,/ (6ir 2 ), 
where is given by Eq. ( fTOb . Then, using the stability crite- 
rion given in Eq. (O, we construct the phase diagrams. For its 
simplicity, we mainly present our results for the equal mass 
(m^ = m^) mixtures, but we also briefly mention the effects 
of mass anisotropy 7^ m; on the phase diagrams. 
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FIG. 1 . (Color online) Saddle-point solutions for the amplitude of the 
order parameter | A | , chemical potential sum /i+ = (/i-f + m ) /2, and 
the chemical potential difference /i_ = (/i^ — /ij.)/2 are shown (in 
units of ef.j.) as a function of the effective s-wave scattering param- 
eter l/(kF,ia e f /). (/i- is shown only for weak interactions where 
it is relevant). This data corresponds to the case where = mj_, 
tf — €f,i and kF.idz — 0.1. The arrow shows the location of 
the topological gapless superfiuid to gapped superfiuid transition dis- 
cussed in the text. 



A. Normal-Superfluid transition 

Using kp.i and cf,i = k 2 F J (27714,) as our length and en- 
ergy scales, respectively, we solve Eqs. ( 1A1I ). ( IA2b and dA3l ) 
numerically. For instance, in Fig. Q] we show self-consistent 
solutions of |A|, /i + and /i_ as a function of l/(k F ^a e ff), 
when = toj., tf = e F ,i and k F ^d z = 0.1. When 
the scattering parameter is smaller than a critical value, i.e. 
V(k-F>4.°e//) —0.32, the saddle-point solution = /i + — 
/j_ — > 1 (in units of e F .]} since |A| = 0, indicating that the 
mixture is a normal Fermi gas i24fl . Beyond this critical value, 
the superfiuid order parameter | A | becomes nonzero indicat- 
ing a quantum phase transition from the normal to a superfiuid 
phase. 
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This transition can be understood from earlier works on 
Cooper pairing with mismatched Fermi surfaces. For in- 
stance, in the case of Fermi gases with unequal densities in 
purely three-dimensions, a superfiuid to normal phase tran- 
sition has been recently observed beyond a critical density 
difference depending on the value of the scattering param- 
eter lfl5l4l8ll . This transition occurs when the difference in 
the chemical potentials (or Fermi momenta) reaches what is 
known as the Clogston-Chandrasekhar limit l25l l26ll . In our 
case, the main mechanism is the same. As can be extracted 
from Eq. ( fTTT i. the mismatch is inevitable in some parts of the 
k-space even for the equal density mixtures considered in this 
manuscript. Therefore, it is energetically more favorable for 
the mixture to be in the normal phase until a critical scattering 
parameter is reached, beyond which Cooper pairing is possi- 
ble. 



B. Topological gapless-gapped superfluidity transition 

With further increase in the scattering parameter, |A| in- 
creases quite rapidly and fi + decreases, with a kink in the 
former quantity at l/(kp^a e ff) ~ —0.22. (We also expect 
a weak kink in /i+ at the same point, but it is not clearly 
seen in the data.) Therefore, the BCS-BEC evolution in 
mixed-dimensional Fermi gases is nonanalytic, i.e. it is not 
a crossover. Recall that, in usual three-dimensional mixtures, 
the evolution of |A| and [i + is analytic for all l/(kp^a e ff), 
and the evolution is just a crossover. The kink in | A| is more 
pronounced for lower values of and it signals a topological 
quantum phase transition as discussed next. 

The excitation spectrum of quasiparticles are determined 
by energies E^i and E^2- At k-space points, the condi- 
tion i?k,s = defines Fermi surfaces of quasiparticles in 
momentum space where the quasiparticle excitation spectrum 
changes from a gapped to a gapless phase. These changes 
in the Fermi surfaces of quasiparticles are topological in na- 
ture ll27ll . and we identify topological quantum phase transi- 
tions associated with the disappearance or appearance of mo- 
mentum space regions of zero quasiparticle energies when ei- 
ther \/(kF^a e ff), if, and/or kp^d z is changed. Note that the 
topological transition occurs without changing the symmetry 
of the order parameter as the Landau classification demands 
for ordinary phase transitions. 

We illustrate the gapless superfiuid phase in Fig. |2j where 
contour maps of E^ s are shown as a function of k z and k y in 
the k x = plane, when m-j- = m±, tf — ep^, kp^d z = 0.1 
and l/(kF^a e ff) = —0.25. In Fig.Q] this data corresponds 
to a point that is slightly on the left hand side of the transi- 
tion point indicated by an arrow. The excitations are gapless 
(£k,s < 0) in the white regions. 

The topological transition could be potentially observed 
through the measurement of the momentum distribution 
of the fermions [28], which can be extracted from Eqs. (|7]i 
and (0. For k-space regions where i?k,i > and £^,2 > 0, 
the corresponding momentum distributions are equal n^ f — 
?ik,i- However, when < and Ek,2 > 0, then rtk.t = 1 
and nk,j. = 0. Similarly, when i?k,i > and E^ 2 < 0, then 
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( b ) E (k x =0,k y ,k z ),2 




-1 1 

k y /k F ,{ 

FIG. 2. (Color online) Contour maps of (a) the quasiparticle excita- 
tion energy -Ek.ii and (b) the negative of the quasihole excitation en- 
ergy -Ek.2 are shown as a function of momentum components k z (in 
units of 1 /dz) and k y (in units of fef ,4.) in the k x — plane. This data 
corresponds to the case where m-f = m;, tf = £f,i, kp^d z — 0.1 
and l/(fcF,4.a e //) = —0.25. The excitations are gapless (i5k,s < 0) 
in the white regions. 

"k,t = an d n Ki = 1- W e illustrate these cases in Fig. [3] 
for the parameters of Fig. [2] Note that, although there are ex- 
cess (or unpaired) j" or I fermions in different regions of the 
k-space, i.e. the bright yellow regions, there are equal number 
of j" and 4. fermions in total. The size of yellow regions may 
look very different in (a) and (b) due partly to the difference 
in scaling factors in k z and k y . 

This topological transition is quantum (T = 0) in nature, 
but its signatures should still be observed at finite tempera- 
tures within the quantum critical region, where the momentum 
distributions are smeared out due to thermal effects. In ad- 
dition, while thermodynamic quantities such as atomic com- 
pressibility, specific heat, and spin susceptibility have power- 
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FIG. 3. (Color online) Contour maps of the momentum distributions 
(a) rik.t for the f fermions, and (b) nk.j. for the 4- fermions are shown 
as a function of momentum components k z (in units of l/d z ) and k y 
(in units of fef,4.) in the k x = plane. This data corresponds to the 
case considered in Fig. [2] where rrif = mi, tf = £f,4,> kp^d z =0.1 
and \/{kF,] r a e ff) = —0.25. Note that when _Ek,i > and £k,2 > 
0, then rik,t = Wk4- However, when i5k,i < and _Bk,2 > 0, then 
rik.t = 1 an d nk 4 — 0> an d wnen -Ek.i > and _Bk,2 < 0, then 
nw,f = and rik, 4. = 1. The densities are 1 (0) in the bright yellow 
(black) regions. 



law dependences on the temperature in the BCS side, they 
have exponential dependences on the temperature and the 
minimum energy of quasiparticle excitations in the BEC side, 
again signaling the existence of a quantum phase transition 
at T = 0. Having discussed the topological classification of 
possible superfluid phases, we are ready to present the saddle- 
point phase diagrams at T = 0, including the stability analysis 
(positive curvature criterion) of the solutions. 



C. Ground state phase diagrams 

In Fig. |U ground state phase diagrams are shown as a 
function of the tunneling amplitude tf and the effective s- 
wave scattering parameter l/(fci?4a e //) for fixed values of 
(a) kF,id z = 0.1, and (b) kF,id z — 1. We indicate nor- 
mal (N), gapped uniform superfluid (SF), gapless uniform 
superfluid (gSF), and phase-separated (PS) regions. The 
normal phase is characterized by a vanishing order parameter 
(A = 0), while the gapped superfluid and gapless superfluid 
phases are both characterized by A > and d 2 flo/dA 2 > 0, 
but with distinct k-space topologies as discussed above. In 
the gSF phase, the unpaired j" and 1 fermions coexist with 
the paired (or superfluid) ones in different k-space regions, 
but there are no unpaired fermions in the SF phase, i.e. all f 
and l fermions are paired. The gSF phase is in some ways 
similar to the Sarma state found in mixtures with unequal den- 
sities fl2ll . but in our case, the gapless superfluid phase is 
unpolarized and most importantly it is stable against phase 
separation. The phase-separated region is characterized by 
d 2 flo/dA 2 < 0, but this region could also be of the FFLO- 
type superfluid having spatial modulations ll29l l3(ill . Such a 
possibility is not considered in this manuscript, and it is left 
as an important problem to address in the future. 

We can understand these phase diagrams as follows. For a 
fixed tf, when the scattering parameter is smaller than a crit- 
ical value, the potential energy is not sufficient to cause pair- 
ing due to mismatch of the Fermi surfaces, and the mixture 
is a normal Fermi gas with A = 0. As shown in Fig. Ufa), 
the critical scattering parameter decreases with increasing tf, 
since increasing tf decreases the mismatch for lower values 
of tf. In the normal region, when tf — 0, Fermi surface 
of the t fermions is a cylindrical shell in the k-space with 
height kpz.f = 27r/<i z and radius kp : f- However, in the low- 
filling limit of I fermions, Fermi surface of the 4- fermions 
is more like a spherical shell with radius kp.i- Note that for 
equal mass and equal density mixtures considered here, the 
k-space volumes enclosed by the cylindrical and spherical 
Fermi surfaces must be equal. Therefore, at tf = 0, there 
is a large mismatch between the two Fermi surfaces when 
kp.f <C fcp,4. -C Tr/d z , and increasing tf from zero decreases 
kp z ,f and increases the ratio kp.f/kp^. When this happens, 
Fermi surface of the j" fermions looks like a prolate spheroid 
(like an american football). This decreases the mismatch for 
small values of tf as long as kp.f < ^fU> an d it 1S quali- 
tatively what happens along the normal-superfluid transition 
boundary in Fig. Ufa), when kp^d z = 0.1 <C 7r. 

In Fig. |4fb), we show the same phase diagram for a higher 
value of kp,id z = 1. Similar to Fig. |4ja), the critical scatter- 
ing parameter for the normal-superfluid transition decreases 
initially with increasing tf, since increasing tf decreases the 
mismatch for lower values of tf. In contrast, beyond a critical 
value of tf [e.g. tf > ep.i in Fig. |4jb)], the critical scattering 
parameter increases with increasing tf. This is because, since 
increasing tf decreases kp. z and increases the ratio kp,f/kp,i, 
it eventually increases the mismatch of the two Fermi surfaces 
after kF,f > kp,i- When this happens, the Fermi surface of 
the j" fermions changes from a prolate spheroid to an oblate 



6 



1.5 



^ 1 



0.5 




(a) kp.xd z = 0.1 


l i|i 
1 1 1 

\ ^11 
\ b£, \ / 

\ *\ / 


/PS 


O 1 


\ *\f 
\ *\ 


SF 


ii ! . 
+ i 


N \ •! 




3. 1 


\ A 

\ j 







2 
1.5 



^ 1 



-0.5 -0.25 0.25 0.5 0.75 

l/(k F ,ia ef f) 

(b) k Rl d z = 1 



0.5 



j4 




i 
i 
i 

O 1 


gSF 


SF 


II \ . 


N 







-2.4 -1.6 -0.8 0.8 

V(^,4 a e//) 

FIG. 4. (Color online) Ground state phase diagrams are shown as 
a function of the tunneling amplitude tf (in units of an d the 
effective s-wave scattering parameter l/(feF,j.o e //). Here, we con- 
sider equal density (n-f = nj,) and equal mass (m-f- = mj,) mixtures, 
for fixed values of kF,i.d z = 0.1 in (a), and kF,id z — 1 in (b). We 
show normal (AT), gapped uniform superfluid (SF), gapless uniform 
superfluid (gSF), and phase-separated (PS) regions. The PS re- 
gion in (b) is close to the normal-superfiuid transition boundary, but 
it is very small and not shown. There is no phase transition across 
fi+ = (green-dashed) line, it is only shown as a reference point. 



spheroid (a disk-shaped ellipsoid). For the case considered in 
Fig.|4ja) where kp,fd z — 0.1, we found that this occurs be- 
yond tf > 100ei?4, but it is not shown. In general, for equal 
mass and equal density mixtures, such a change is expected to 
occur beyond tf > l/(2m,fd 2 z ), and this expectation is con- 
sistent with our numerical findings. 

In addition, in Fig. [4] the solid-black lines correspond to 
the transition boundary between the gapless superfluid (gSF) 
and gapped superfluid (SF) phases. However, the stability 
criterion given in Eq. (0 is not satisfied in some parts of the 
gapless superfluid region, indicating a phase separation. This 
occurs between the dashed-blue and the solid-black lines in 
Fig.Ua). In contrast, phase separation occurs in a tiny region 
very close to the normal-superfiuid transition boundary (solid- 
red line) in Fig.|4|b), and it is not shown. The gapless super- 
fluid phase phase is in some ways similar to the Sarma state 



found in mixtures with unequal densities 1 12], but in our case, 
gSF is unpolarized and most importantly it is stable against 
phase separation in a considerably large region as shown in 
the figures. In this way, our gSF phase is similar to those 
of Refs. lfl3l [l4ll . which are recently proposed for ultracold 
atomic systems in other contexts. 

Before concluding, we would like to comment on the phase 
diagram of mixed-dimensional two-species Fermi-Fermi mix- 
tures. When t and J, fermions have different masses, the 
phase boundaries shift left (right) when the j" species is 
heavier (lighter) than the I fermions. For instance, in the 
case of 6 Li- 40 K mixtures, the tf — >• limit of the normal- 
superfiuid boundary for the case considered in Fig. Ufa) shifts 
to l/(kp^a e ff) « —0.6 when rrif = 6.64raj, (when 40 K 
atoms correspond to f) and to 1 / (kF,ia>ef f) ~ 0.6 when 
nif = 0.15m^ (when 6 Li-atoms correspond to f). 



IV. CONCLUSIONS 

Motivated by a very recent experiment involving mixed- 
dimensional Bose-Bose mixtures here we investigated 
the ground state phase diagram of superfluidity for mixed- 
dimensional Fermi-Fermi mixtues in the BCS-BEC evolution. 
In this recent experiment, a species-selective one-dimensional 
optical lattice is applied to a two-species mixture of bosonic 
atoms, such that only one of the species feel the lattice po- 
tential, and is confined to a quasi-two-dimensional geometry, 
while having negligible effect on the other, that is leaving it 
three dimensional. We considered a similar problem with two- 
species mixtures of fermionic atoms, where both species are 
confined to quasi-two-dimensional geometries determined by 
their hoppings along the lattice direction. 

We considered equal-density mixtures at zero temperature, 
and after solving the saddle-point self-consistency equations, 
we constructed the phase diagrams using some stability crite- 
rion. We found normal, gapped superfluid, gapless superfluid, 
and phase separated regions. The gapped superfluid and gap- 
less superfluid phases are identified with the disappearance or 
appearance of momentum space regions of zero quasiparticle 
energies. In particular, we found a stable gapless superfluid 
phase where the unpaired j" and I fermions coexist with the 
paired (or superfluid) ones in different momentum space re- 
gions. This phase is in some ways similar to the Sarma state 
found in mixtures with unequal densities [12], but in our case, 
the gapless superfluid phase is unpolarized and most impor- 
tantly it is stable against phase separation. We also argued 
that the topological transition from the gapped superfluid to 
the gapless superfluid could be potentially observed through 
the measurement of the momentum distribution l28ll . 

There are several ways to extend this work. First, the possi- 
bility of FFLO-type superfluid phases [29, 30], where Cooper 
pairs have finite center of mass momentum, leading to a spa- 
tially modulated phase, is not considered in this manuscript. 
This is an important problem to address due to its relevance 
to condensed-matter systems. Second, our calculation is 
based on the saddle-point self-consistency equations, which 
are known to be sufficient to qualitatively describe the entire 
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BCS-BEC evolution, at least for the usual three-dimensional 
mixtures at low temperatures. However, corrections beyond 
the saddle point could be important in stabilizing or destabliz- 
ing especially the gapless superfluid phase [31]. Third, we 
used a single-band model to describe the optical lattice po- 
tential, and the effects of higher bands could become impor- 
tant near the strongly interacting regime. Lastly, atomic sys- 
tems are not uniform since confining trapping potentials are 
always present, and finite-size effects due to such potentials 
could also be analyzed. 
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Appendix A: Self-consistency equations at T = 

For numerical purposes, the self-consistency equations can 
be simplified as follows. At zero temperature, since the Fermi 
function f(x) turns into a heaviside step function 9(—x), 
Eqs. ©, (0 and ^ become 



m + V 



Ana 



eff 



E 



1 



1 - 



-E k ,i) 



-£k, 2 ) 



2e k , 



2E U 



N t = [K| 2 #(-£k,i) + \v k \ 2 6(E K2 )] , (A2) 

k 

N l = E D w k| 2 #(-£k, 2 ) + k| 2 #(£k,i)] , (A3) 



where the k-space sums are, ^ k = [V/ (2ir) 3 } J <i 3 k = 

[V^/(27r 2 )] Jq dk z J °° k±dk±. Note that pairing occurs 
only in the k-space regions where both Ek i and E^2 nave 
the same (positive) sign. When _Ek,i < and £^,2 > or 
vice versa, the first term (on the right hand side) of Eq. ( lAlb 
inside the parentheses vanishes, reflecting that the pairing is 
not allowed for those k-space regions, and the quasiparticle 
and quasihole excitations are gapless. 

In order to perform the integration over k± by hand, we 
need to find the k-space regions where the excitations are gap- 
less, i.e. E'k.i < or £^,2 < 0. The zeros of E'k.i and £k,2 
are determined by real and positive solutions of 







m + /m_ 



)x 2 + 2(£ fcai+ - m+£k„,-/m-)x 



+ £fc*,t£fc*,; + |A| 



(A4) 



where x — fc+/(2m+) and £k x ,± = e k z ,± — M±- Here, 
we introduced e kz ,± = (ek m ,f ± efc 2 ,j.)/2, where e k ^ a = 
2t a [1 — cos(k z d z )) is the energy dispersions in the k z direc- 
tion. Solutions of this equation (a;< and x>) depend on k z , 
and they give the locations of the zeros of £^,1 and £k.2 
in the fc^-axis as a function of k z . For instance, if both 
x< and cc> are real and positive, then Ek.i is gapless for 
x < < k\/(2m + ) < x> in some k z region z< t i < k z < z>,i, 
and i?k,2 is gapless for x< <k\j (2m+) < x > in some other 
k z region z<,2 < k z < 2>,2- If only x> is real and positive, 
£k,i is gapless for < k\/(2m + ) < £> in some k z region 
z<.i < k z < z>.i, and i?k,2 is gapless for < k\/{2m + ) < 
Xy in some other k z region z< 2 < k z < z>2- If there is 
no real and positive solution, then the excitations are always 
gapped. 

Given the gapless k-space regions, Eq. ( IAU can be written 

as 



2 eff 



z >,i r z >* 2 



Z <,1 z < ,2 



dk z In 



x> + Ck.,+ + y/{x>+Ck„+) 2 + \A\ 



Txjd z 



dk z In 



2e kz 



where the first term on the right hand side is coming from the gapless, but the second term is from the gapped k-space regions. 
Similarly, Eq. (lA2t can be written as 

47r 2 n-j. 



771+ 



(/■ • 


r z >,2 


J dk z 


x> 




J z <,2 







Z >,1 



x> - \ (x> +6,,+) 2 + |A| 2 -x< + J( X< +6„+) 2 + |A 



TT/d z 



2 1 dk z (x > -x < )+ I dfc, + |A| 3 -&,,+ ), 



(A6) 



where again the first two terms on the right hand side is com- 
ing from the gapless, but the third term is from the gapped k- 
space regions. The density of J, fermions can be obtained by 
substituting 1 — > 2 in the integration limits of the second term. 
Our numerical calculations show that integrating k± by hand 
as described above and calculating the remaining k z integral 
numerically is a much more stable approach compared to the 



one where both integrations are calculated numerically. In 
particular, this approach converges much faster than the latter 
on the BCS side, where integrations involve gapless k-space 
regions. 
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Appendix B: Normal-Superfluid Phase Boundary at T = 



The phase boundary for the normal-superfiuid transition 
can be found from Eqs. ([6j, (0 and ([H} by setting A = 0. 
Therefore, at the transition boundary, the self-consistency 
equations are uncoupled, i.e. the gap equation determines the 
critical effective scattering length, and the chemical potentials 
are determined by the number equations. At zero temperature, 



this leads to 

771+ V 



y 



i 



i - 0(-&, t ) - 0(-6a) 



2ek,+ 



(Bl) 
(B2) 



Note that Eq. ( IB2b is the number equation for noninteracting 
fermions, and it is already solved in Sec. Ill CI Similar to Ap- 
pendix A, we can also simplify Eq. © by finding the zeros of 
£k,er- This leads to 



- — = dk z In f ^ ) + In ( £fc -+ ) + r 7 "" dk z In (f**±) , (B3) 



where fc Fz!< = min{/c F2iT , fc Fz ,;} and /c FZj> = 
max{fc FZ; -f-, fc FZl j,}. For instance, when fc Fz ,< = kp z ,^ then 
<=4- and >=t in the integrands. Here, we again emphasize 
that integrating k± by hand as described above, and calcu- 



lating the remaining k z integral numerically is a much more 
stable approach compared to the one where both integrations 
are calculated numerically. 
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